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Q ' Abstract 

(N ! 

We present the perturbative parts of the structure functions F% and F£ for a gluon 
target having nonzero transverse momentum squared at order a s . The results of the 
<^ ■ double convolution (with respect to the Bjorken variable x and the transverse momen- 

tum) of the perturbative part and the unintegrated gluon densities are compared with 
| HERA experimental data for F%. The contribution from F£ structure function ranges 

10 ~r~ 30% of that of F£ at the HERA kinematical range. 
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1 Introduction 
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Recently there have been important new data on the charm structure function (SF) F 2 C , 
of the proton from the HI |l| and ZEUS Collaborations at HERA, which have probed 
the small-x region down to x = 8 x 10~ 4 and x = 2 x 10~ 4 , respectively. At these values of 
x, the charm contribution to the total proton SF, F^, is found to be around 25%, which is a 
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considerably larger fraction than that found by the European Muon Collaboration at CERN 
PI at larger x, where it was only ~ 1% of F%. Extensive theoretical analyses in recent years 
have generally served to confirm that the _F 2 C data can be described through perturbative 
generation of charm within QCD (see, for example, the review in Ref. [Q and references 
^ . therein). 

We note, that perhaps more relevant analyses of the HERA data, where the x values are 
quite small, are those based on BFKL dynamics ||, because the leading ln(l/x) contributions 
are summed. The basic dynamical quantity in BFKL approach is the unintegrated gluon 
distribution (p g (x, kjj [] (f g is the (integrated) gluon distribution multiplied by x and k± is 
the transverse momentum) 



f g (x,Q 2 ) = f g (x,Ql)+ /jj dk 2 ± <p g (x,k 2 ± ) 



o ~ ' " " ' 



which satisfies the BFKL equation. The integral is divergent at the lower limit and it leads 
to the necessity to use the difference f g (x, Q 2 ) — f g (x, Qq) with some nonzero Q^. 

Then, in the BFKL-like approach (hereafter the fc r factorization approach || [7J is used) 
the SFs F 2 c L (x, Q 2 ) are driven at small x by gluons and are related in the following way to 
the unintegrated distribution (p g (x, k±): 

Fl L (x,Q 2 ) = f- I dk 2 ± Cl L (z,Q 2 ,mlk 2 ± ) <p g {x/z,kti, (2) 

J x Z J 



hereafter and are the hadron and the gluon 4-momentums, respectively, and is the photon 
4-momentum. 
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The functions Cf L (x, Q 2 , m 2 ., fc 2 ) may be regarded as the structure functions of the off- 
shell gluons with virtuality k\ (hereafter we call them as hard structure functions 2 ) . They 
are described by the quark box (and crossed box) diagram contribution to the photon-gluon 
interaction (see Fig. 1 in ||). 

The purpose of the paper is to present the results for these hard SF C% L (x, Q 2 , m 2 ., kj_) 
and to analyze experimental data for F%(x, Q 2 ) by applying Eq. (p|) with different sets of 
unintegrated gluon densities (see Ref. 0, |TU|) and to give predictions for the longitudinal 
SF F[(xb, Q 2 )- 

It is instructive to note that the results should be similar to those of the photon-photon 
scattering process. The corresponding QED contributions have been calculated many years 
ago in Ref. 0] (see also the beautiful review in Ref. Hl2|). Our results have been calculated 
independently in || (based on approaches of |13|, [L4[]) and they are in full agreement with 



TTJ. However, we hope that our formulas which are given in a simpler form could be useful 
for others. 



2 Hard structure functions 



The gluon polarization tensor (hereafter the indices a and (3 are connected with gluons), 
which gives the main contribution at high energy limit, has the form: 
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where /3 2 = 1 - 46a; 2 , b = -k 2 /Q 2 = k 2 JQ 2 > 0. 

Contracting the corresponding photon projectors, we have: 
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where the normalization factor /C = e 2 a s {Q 2 ) / '(47r) x, a — m 2 /Q 2 , e c = 2/3 is charm quark 
charge and 
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1 - 1 - 2x{l + b - 2a) [1 - x(l + b + 2a)] /i 
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+ (2o-6)(l -2a)x 2 / 2 
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;i - (1 + b)x) - 2x(bx(l - (1 + b)x)(l + b-2a) + a/3 2 ) / x 



+ bx 2 (l- (l + b)x)(2a -b) f 2 
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2 x(l-x(l + 6))(l + 6-2a) + a/3 2 / 



2 This notation reflects the fact that structure functions F£ L connect with the functions C| L at the same 



form as cross-sections connect with hard ones (see [pi W). 
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- x(l-x(l + 6))(l-2a) h 
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2-(l + 26x 2 ) h-bx 2 f 2 
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(l-(l + 6)x)' 
For the important regimes when k 2 
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with ones of Refs. [0, [TJ|. Notice that our results in Eq. 
in Ref. 



1 - /? 2 ^ 2 

0, the results are coincided 
should be also agree with those 



16| but the direct comparison is quite difficult because the structure of their results 



is quite cumbersome (see Appendix A in Ref. fll6|). We have found numerical agreement in 
the case of F 2 (x,Q 2 ) for several types of unintegrated gluon distributions (see Fig. 4 in M). 



3 Comparison with F 2 C experimental data and predic- 



tions for F? 



With the help of the results obtained in the previous Section we have analyzed HERA 
data for SF F% from ZEUS [[| collaboration. 

Notice that in Ref. the k\ -integral in the r.h.s. of Eq.(0) has been evaluated using the 
BFKL results for the Mellin transform of the unintegrated gluon distribution and the Wilson 
coefficient functions have been calculated for the full perturbative series at asymptotically 
small x values. Since we would like to analyze experimental data for the SF F 2 C , we have an 
interest to obtain results at quite broad range of small x values. For the reason we need in 
a parameterization of unintegrated gluon distribution function. 

We consider two different parametrizations for the unintegrated gluon distribution (see 



]): the Ryskin-Shabelski (RS) one |17] and the Blumlein one [|TJ 



In Fig. 1 we show the SF F% as a function x for different values of Q in comparison 
with ZEUS § experimental data. We see that at large Q 2 (Q 2 > 10 GeV 2 ) the SF F 2 C 
obtained in the factorization approach is higher than the SF obtained in the standard 
parton model with the GRV [I9| gluon density at the LO approximation (see curve 1) and 



has a more rapid growth in comparison with the standard parton model results, especially 
at Q 2 ~ 130 GeV 2 [0. At Q 2 < 10 GeV 2 the predictions from perturbative QCD (in GRV 
approach) and those based on the hp factorization approach are very similar [] and show 
a disagreement with data below Q 2 = 7 GeV 2 0. Unfortunately the available experimental 
data do not permit yet to distinguish the k? factorization effects from those due to boundary 
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The results for the SF F£ obtained in perturbative QCD and from the kj> factorization 
approach are quite similar to the F£ case discussed above. The ratio R c = F[/F£ is shown 
in Fig. 2. We see that R c ~ 0.1 -j- 0.3 in a wide region of Q 2 . The estimation of R c is very 

3 This fact is due to the quite large value of Q§ = 4 GeV 2 chosen here. 

4 A similar disagreement with data at Q 2 < 2 GeV 2 has been observed for the complete structure function 
F 2 (see, for example, the discussion in Ref. ElJ and reference therein). We note that the insertion of higher- 
twist corrections in the framework of usual perturbative QCD improves the agreement with data (see Ref. 
p2]) at quite low values of Q 2 . 
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Figure 1: The structure function F£(x, Q 2 ) as a function of x for different values of Q 2 
compared to ZEUS data [§. Curves 1, 2 and 3 correspond to SF obtained in the standard 
parton model with the GRV Jl9| gluon density at the leading order approximation and to SF 
obtained in the kr factorization approach with RS |l7j and Blumlein (at Qq = 4 GeV 2 ) fl8H 
parameterizations of unintegrated gluon distribution. 
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close to the results for R = Fl/(F 2 — Fl) ratio (see Refs. [p3fl-[p5|). We would like to note 



that these values of R c contradict the estimation obtained in Refs. |2|, [j]. The effect of R c 
on the corresponding differential cross-section should be considered in the extraction of _F 2 C 
from future more precise measurements. 

For the ratio R c we found quite flat x-behavior at low x in the low Q 2 region (see 
Fig. 2), where approaches based on perturbative QCD and on kx factorization give similar 
predictions (see Fig 1). It is in agreement with the corresponding behaviour of the ratio 
R = F L /(F 2 - F L ) (see Ref. @]) at quite large values of A P [] (A P > 0.2 - 0.3). The low x 
rise of R c at high Q 2 disagrees with early calculations [2~3| in the framework of perturbative 



QCD. It could be due to the small x resummation, which is important at high Q 2 (see Fig 
1). We plan to study this effect in future. 

4 Conclusions 

We have presented the results for the perturbative parts of the SFs F| and F£ for a gluon 
target having nonzero momentum squared, in the process of photon-gluon fusion. 

We have applied the results in the framework of k? factorization approach to the analysis 
of present data for the charm contribution to F 2 and we have given the predictions for 
F[. The analysis has been performed with two parameterizations of unintegrated gluon 
distributions for comparison. We have found good agreement of our results with experimental 
HERA data for F 2 C , except at low Q 2 (Q 2 < 7 GeV 2 ) Q. We have also obtained quite large 
contribution of the SF F£ at low x and high Q 2 (Q 2 > 30 GeV 2 ). 

We would like to note the good agreement between our results for F% and the ones 
obtained in Ref. |29|] by Monte-Carlo studies. Moreover, we have also good agreement with 
fits of HI and ZEUS data for F 2 (see recent reviews in Ref. [[[IJ and references therein) 
based on perturbative QCD calculations. But unlike to these fits, our analysis uses universal 
unintegrated gluon distribution, which gives in the simplest way the main contribution to 
the cross-section in the high-energy limit. 

It could be also very useful to evaluate the complete F 2 itself and the derivatives of F 2 
with respect to the logarithms of 1/x and Q 2 with our expressions using the unintegrated 
gluons. 

The consideration of the SF F 2 in the framework of the leading-twist approximation of 
perturbative QCD (i.e. for "pure" perturbative QCD) leads to very good agreement (see Ref. 
| |2~Tf and references therein) with HERA data at low x and Q 2 > 1.5 GeV 2 . The agreement 
improves at lower Q 2 when higher twist terms are taken into account [22]. As it has been 
studied in Refs. |2"I| , ^2]] , the SF F 2 at low Q 2 is sensitive to the small- £ behavior of quark 
distributions. Thus, our future analysis of F 2 in a broader Q 2 range in the framework of 
kj- factorization should require the incorporation of parametrizations for unintegrated quark 
densities, introduced recently (see Ref. |JT) and references therein). 

The SF Fl is very interested too (see J32|). The structure function Fl depends strongly 
on the gluon distribution (see, for example, Ref. |33|, 0), which in turn is determined [|35]] by 
the derivative dF 2 /dlnQ 2 . Thus, in the framework of perturbative QCD at low x the relation 

5 The behaviour is in agreement with previous studies ^3). Note that at small values of Ap, i.e. when 
x~ Ap ~ Const, the ratio R has got strong negative NLO corrections (see Refs. [^7] ^8|) and tends to zero 
at x — > after some resummation done in p8f . 

6 It must be noted that the cross section of inelastic cc— and bb— pair photoproduction at HERA are 
described by the Blumlein parametrization at a smaller value of Qq (Qq = 1 GeV 2 ) (see 0). 
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between Fl, F2 and (ii^/rflnQ 2 could be violated by non-perturbative contributions, which 
were expected to be important in the F L case (see Refs. [^6], |37j). The application in |32 



of present analysis to F L gives a "non pure" perturbative QCD predictions for the structure 
function which are in good agreement with data p5|| and with the "pure" perturbative results 



of Ref. p3|, that determines a small value of nonperturbative corrections to SF Fl at low x 



and quite large Q 2 values. 

Acknowledgments. Authors would like to express their sincerely thanks to the Orga- 
nizing Committee of the XVIth International Workshop "High Energy Physics and Quantum 
Field Theory" for the kind invitation, the financial support at such remarkable Conferences, 
and for fruitful discussions. 

A.V.K. was supported in part by Alexander von Humboldt fellowship and INTAS grant 
N366. CP. acknowledges the support of Galician research funds (PGIDTOO PX20615PR) 
and Spanish CICYT (FPA2002-01161). N.P.Z. also acknowledge the support of Royal 
Swedish Academy of Sciences. 

References 



[i] 

[2] 
[3 

[4] 
[5] 



[6] 
[7] 



[9 
[10 

[11 
[12 

[13 
[14 

[15 

[16 
[17 



HI Coll.: S. Aid et al, Z.Phys. C72 (1996) 593; Nucl.Phys. B545 (1999) 21. 

ZEUS Coll.: J. Breitweg et al, Phys.Lett. B407 (1997)402; Eur.Phys.J. C12 (2000)35. 

EM Coll.: J.J. Aubert et al, Nucl.Phys. B213 (1983) 31; Phys.Lett. B94 (1980) 96; 
B110 (1983) 72. 

A.M. Cooper-Sarkar et al, Int.J.Mod.Phys. A13 (1998) 3385. 

L.N. Lipatov, Sov. J.Nucl.Phys. 23 (1976) 338; E.A. Kuraev et al, Sov.Phys.JETP 44 
(1976) 443, 45 (1977) 199; Ya.Ya. Balitzki, L.N. Lipatov, Sov. J.Nucl.Phys. 28 (1978) 
822; L.N. Lipatov, JETP 63 (1986) 904. 

S. Catani et al, Phys.Lett. B242 (1990) 97; Nucl.Phys. B366 (1991) 135. 

J.C. Collins, R.K. Ellis, Nucl.Phys. B360 (1991) 3; E.M. Levin et al, Sov. J. Nucl. Phys. 
53 (1991) 657. 

A.V. Kotikov et al, Preprint US-FT/7-01 Qhep-ph/01 071351) . 



A. V. Lipatov et al., Mod.Phys.Lett. A15 (2000) 695; 1727. 



Bo Andersson et al., [hep-ph/0204115 



V.N. Baier et al., JETP 50 (1966) 156; V.G. Zima, Yad.Fiz. 16 (1972) 1051. 
V.M. Budnev et al., Phys.Rept. 15 (1975) 181. 

D. I. Kazakov, A.V. Kotikov, Theor.Math.Phys. 73 (1987) 1264; Nucl.Phys. B307 (1988) 
721; Nucl.Phys. B345 (1990) 299. 

A.V. Kotikov, Theor.Math.Phys. 78 (1989) 134; Phys.Lett. B375 (1996) 240; Phys.Lett. 
B254 (1991) 158; Phys.Lett. B259 (1991) 314; Phys.Lett. B267 (1991) 123; in: Proc. of 
the XXXV Winter School, Repino, S'Peterburg, 2001 flhep-ph/0112347|) . 

E. Witten, Nucl.Phys. B104 (1976) 445; M. Gluck, E. Reya, Phys.Lett. B83 (1979) 98; 

F. M. Steffens et al, Eur.Phys.J. Cll (1999) 673. 

G. Bottazzi et al, JHEP 9812 (1998) 011. 

M.G. Ryskin, Yu.M. Shabelski, Z.Phys. C61 (1994) 517; C66 (1995) 151. 



7 



[18] J. Blumlein, Preprints DESY 95-121 (|hep-ph/9506403|) ; DESY 95-125 flKep"- 
| ph/9506446]) . 

[19] M. Gluck et al., Z. Phys. C67 (1995) 433. 

[20] A.V. Lipatov, N.P Zotov, in Proc. of the 8th Int. Workshop DIS 2000 (2000), World 
Scientific, p. 157. 

[21] A.V. Kotikov, G. Parente, Nucl.Phys. B549 (1999) 242; |hep-ph/001035l ; in Proc. of 
the Int. Conference PQFT98 (1998), Dubna (|hep-ph/ 98 10223; ); in Proc. of the 8th 
Int. Workshop on Deep Inelastic Scattering, DIS 2000 (2000), Liverpool, p. 198 ( |hep- 

| ph/0006T97D ; Preprint US-FT/3-02 ( |hep-ph/ 0207276]) . 

[22] A.V. Kotikov, G. Parente, in Proc. Int. Seminar Relativistic Nuclear Physics and Quan- 
tum Chromodynamics (2000), Dubna ( [hep-ph/00 12299]) ; in Proc. of the 9th Int. Work- 
shop on Deep Inelastic Scattering, DIS 2001 (2001), Bologna (|hep-ph/ 01 06175) ). 

[23] A.V. Kotikov, JETP 80 (1995) 979; A.V. Kotikov, G. Parente, in Proc. Int. Workshop 
on Deep Inelastic Scattering (1996), Rome, p. 237 flhep-ph/9608409|) ; Mod.Phys.Lett. 
A12 (1997) 963; JETP 85 (1997) 17; |hep-ph/9609239| . 

[24] R.S. Thorne, Phys.Lett. B418 (1998) 371. 

[25] HI Coll.: S. Aid et al, Phys.Lett. B393 (1997) 452; CCFR/NuTeV Coll.: U.K. Yang et 
al., |hep-ex/001000H A. Bodek, |hep-ex/001050E| 7 



[26] A.V. Kotikov et al., Theor. Math. Phys. 84 (1990) 744; Theor. Math. Phys. Ill (1997) 
442; A.V. Kotikov, Phys. Atom. Nucl. 56 (1993) 1276; L. L. Jenkovszky et al., Sov. J. 
Nucl. Phys. 55 (1992) 1224;. JETP Lett. 58 (1993) 163; Phys. Lett. B 314(1993) 421. 

[27] S. Keller et al., Phys.Lett. B270 (1990) 61; L.H. Orr, W.J. Stirling, Phys.Rev.Lett. B66 
(1991) 1673; E. Berger, R. Meng, Phys.Lett. B304 (1993) 318. 

[28] A.V. Kotikov, JETP Lett. 59 (1994) 1; Phys.Lett. B338 (1994) 349. 

[29] H. Jung, Nucl.Phys. (Proc.Suppl.) 79 (1999) 429. 

[30] G. Wolf, Preprint DESY 01-058 (|hep-ex/0105055|) . 

[31] M.A. Kimber et al., Phys.Rev. D63 (2001) 114027. 

[32] S. Catani and F. Hautmann, Nucl.Phys. B427 (1994) 475; A.V. Kotikov et al., |Eepj 

| ph/0207226] . 

[33] A.M. Cooper-Sarkar et al, Z.Phys. C39 (1988) 281; A.V. Kotikov, G. Parente, in: Proc. 
of Madrid Workshop on Low x Physics (1997) pp. 71-76 ( |hep-ph / 9710252] ) . 

[34] A.V. Kotikov, Phys.Atom.Nucl. 57 (1994) 133; Phys.Rev. D49 (1994) 5746. 

[35] K. Prytz, Phys.Lett. B311 (1993) 286; A.V. Kotikov, JETP Lett. 59 (1994) 667; A.V. 
Kotikov, G. Parente, Phys.Lett. B379 (1996) 195. 

[36] J. Bartels et al., Eur.Phys.J. C17 (2000) 121; 

[37] A.V. Kotikov at al., Z. Phys. C58 (1993) 465; G. Parente et al., Phys.Lett. B333 (1994) 
190; A.V. Kotikov, V.G. Krivokhijine, Dubna preprint E2-2001-190 ( |hep-ph/01 082241) ; 
in: Proc. of the XVIth International Workshop "High Energy Physics and Quantum 
Field Theory 11 (2001), Moscow; ( |hep-ph/020622lD ; in: Proc. of the Int. Workshop 
"Renormalizatwn Group 2002" (2002), High Tatras, Slovakia ( |hep-ph/ 0207222]) . in: 
Proc. of the X Int. Workshop on Deep Inelastic Scattering (DIS2002) (2002), Cracow 
( Ihep-ph/OflM&jQ . 



8 



